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The null result of a search for pulsational variations of the surface 
magnetic field in the roAp star 7 Equulei 
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ABSTRACT 

We describe an analysis of the time-resolved measurements of the surface magnetic field in 
the roAp star 7 Equ. We have obtained a high-resolution and high S/N spectroscopic time- 
series, and the magnetic field was determined using Zeeman resolved profiles of the Fe II 
6149.25 A and Fe 1 6173.34 A lines. Contrary to recent reports we do not find any evidence of 
magnetic variability with pulsation phase, and derive an upper limit of 5-10 G for pulsational 
modulation of the surface magnetic field in 7 Equ. 
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1 INTRODUCTION 

After discovery of the conspicuous radial velocity (RV) pulsa- 
tional variations in a sample of rapidly oscillating magnetic pecu- 
liar (roAp) stars (Kanaan & Hatzes 1998 Savanov, Malanushenko 

6 Ryabchikova 1999 Kochukhov & Ryabchikova 1200 la for 

7 Equ; Baldry et al. .l998. Baldry & Bedding 2000l Kochukhov & 



Ryabchikova 12001b| for a Cir and HD 83368), attempts to search 
for magnetic field variations over the pulsational period have been 
made. First, Hubrig et al. ( 2004 ( tried to measure pulsational vari- 
ability of the longitudinal magnetic field in six roAp stars. Their 
sample included 7 Equ ~ probably one of the most favourable stars 
for this kind of investigation. 7 Equ is a bright northern roAp star 
with a strong magnetic field, and with one of the largest pulsational 
RV amplitudes which exceeds 1000 ms~^ in individual spectral 
lines. The extremely slow rotation of 7 Equ, leading to very sharp 
spectral lines, makes this star the best candidate for any study of the 
pulsational variability in spectroscopy. Hubrig et al. <2004> used 
low-resolution Zeeman time-series observations and measured Bi 
using hydrogen lines and unresolved blends of metal lines. They 
failed to detect any variability beyond the formal errors of their 
measurements which were 40-100 G. 

According to a coarse theoretical estimate made by Hubrig 
et al. <2004t there should exist a linear relation between magnetic 
field variability over the pulsational cycle and the RV amplitudes. 
In roAp stars the largest RV amplitudes are observed in the lines of 
first and second ions of rare-earth elements (REE), while they are 
usually below measurement errors for the lines of iron group ele- 
ments. Thus, high-resolution spectroscopy and spectropolarimetry 
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of selected REE spectral lines is a more promising tool for an in- 
vestigation of possible rapid magnetic oscillations in roAp stars. 
Taking this into account Leone & Kurtz (2003^1 obtained a high- 
resolution (R=115 000), high S/N time-series of observations of 
7 Equ with a circular polarization analyzer, and measured Bi using 
four Nd III lines. They reported the discovery of pulsational varia- 
tions with amplitudes between 1 12 to 240 G and, more surprisingly, 
with discrepant phases of magnetic maximum for different Nd III 
lines. Leone & Kurtz's result was based on only 18 time-resolved 
spectra. A year later Kochukhov, Ryabchikova & Piskunov <2004> 
obtained a time-series of polarimetric observations of 7 Equ with 
a smaller resolving power (R=38 000), but acquired more than 200 
spectra over 3 nights, more than compensating for lesser quality 
of individual spectra. Kochukhov et al. <2004> used simultaneously 
13 Ndlll lines for magnetic measurements which allowed them to 
achieve a substantial reduction of the error of the Bi determina- 
tions. They did not confirm longitudinal field variability over the 
pulsational period in 7 Equ and gave a conservative upper limit of 
«40 G for the amplitude of pulsational modulation of Bi deter- 
mined from Nd III lines. 

Another attempt to search for possible rapid magnetic vari- 
ability in 7 Equ was made by Savanov, Musaev & Bondar 1 200^. 
They measured the surface magnetic field Bs variations over the 
pulsational period using the Fell 6149.25 A line observed in un- 
polarized light. Due to a very simple Zeeman pattern (two equally 
separated n- and a-components) this line is ideal for Ba measure- 
ments (see Mathys et al. |1997t . Savanov et al. reported a 1.8(t de- 
tection of Bs variability with an amplitude of 99±53 G. At the 
same time they did not find periodic variations of RV measured for 
the individual Zeeman resolved components of the Fe II 6149.25 A 
doublet exceeding their error limit (100-120 ms~^). The authors 
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Figure 1. Comparison of the average CFHT spectra of 7 Equ (histogram) in the 6149 (left panel) and 6173 A (right panel) regions with magnetic spectrum 
synthesis (sohd hne). The length of the line underlining each identification is proportional to the strength of con'esponding spectral feature. The La II 6172.72 A 
and Eu II 6173.03 Alines were computed taking into account hyperfine splitting. For clarity we do not show individual hfs components, but indicate the range of 
their central wavelengths with horizontal bars. The shaded areas illustrate our selection of the wavelength intervals for the centroid line position measurements. 



Table 1. Journal of spectroscopic observations of 7 Equ. 



UTdate 


Instrument/ 


Wavelength 


Resolution 


Exposure 


Number of 


Start HJD 


End HJD 


Typical 




telescope 


region (A) 


(A/AA) 


time (s) 


exposures 


(2450000+) 


(2450000+) 


SNR 


22/07/1999 


CES/ESO 3.6-m 


6140-6165 


166000 


60 


31 


1381.78344 


1381.82390 


190 


26/09/2002 


Gecko/CFHT 


6104-6194 


115 000 


90 


64 


2543.82314 


2543.92191 


230 



used high-resolution R=120000 time-series observations, but the 
S/N of a single spectrum did not exceed 40-60. Clearly, the result 
of Savanov et al. <2003> is marginal and needs to be confirmed or 
rejected with data of better quality. 

In this paper we present the results of a new search for pulsa- 
tional variations of Bs in 7 Equ using high-resolution and high S/N 
time-resolved observations of this star. These observational data al- 
lowed us to obtain precise measurements of the magnetic field in 
7 Equ, and strongly constrain possible changes of Bs over the pul- 
sation cycle of the star 



2 OBSERVATIONS AND DATA REDUCTION 

The time-resolved observations of 7 Equ were obtained using the 
single-order //4 Gecko coude spectrograph with the EE VI CCD 
at the 3.6-m Canada-France-Hawaii telescope. TableQgives details 
of this spectroscopic time-series dataset. The spectra cover approx- 
imately the spectral window 6104-6194 A. This wavelength inter- 
val contains two Zeeman resolved lines. Fell 6149.25 A and Fei 
6173.34 A. It also has strong Nd III and Prill lines, optimal for in- 
vestigation of pulsational variability of roAp stars, as well as strong 
and weak lines of other elements such as Si, Ca, Cr, and Ba. 

The spectra were reduced using standard IRAF tasks. Each 
stellar, flat and calibration frame had a mean bias subtracted and 
was then cleaned of cosmic ray hits and extracted to one dimension. 
Extracted stellar spectra were divided by an extracted mean flat 
field, and the continuum was fit with a third-order Legendre poly- 
nomial, using the same rejection parameters for all spectra so that 
the continuum fit is as uniform as possible. The wavelength scale 
was established using about 40 lines of a ThAr emission lamp, re- 



sulting in an RMS scatter about the adopted pixel- wavelength poly- 
nomial (a sixth-order Legendre polynomial) of about 5 X 10^* A. 
The wavelength scale was linearly interpolated between ThAr lamp 
spectra taken before and after the stellar series, but the spectra were 
not resampled to a linear wavelength spacing. 

In this paper we also used 31 very high resolution time- 
resolved observations of 7 Equ analysed by Kochukhov & 
Ryabchikova (2001a I. These time-resolved data were obtained 
with the Coude Echelle Spectrograph (CES), fiber-linked to the 
Cassegrain focus of the ESQ 3.6-m telescope. The highest reso- 
lution CES image slicer and the ESQ CCD#38 were used, allowing 
us to reach a resolving power of A/ A A = 166 000 and record spec- 
tra in the 6140-6165 A wavelength interval. We refer the reader to 
Kochukhov & Ryabchikova 1 2001a 1 for other details of the acqui- 
sition and reduction of the CES spectra of 7 Equ. 



3 MAGNETIC FIELD AND RADIAL VELOCITY 
MEASUREMENTS 

As a first step in our analysis of 7 Equ we computed synthetic spec- 
tra in the region of the Zeeman resolved iron lines Fe II 6149.25 A 
and Fel 6173.34 A (hereafter we refer to these spectral features 
simply as the Fell and Fel lines) using abundances, Fe stratifi- 
cation and model atmosphere from the study of Ryabchikova et 
al. (20021. Synthetic spectra of 7 Equ were calculated with the 
SYNTHMAG code (Piskunovri999>, modified to take into account 
vertical stratification of chemical abundances. The simplest model, 
with a constant magnetic field over the stellar surface, was adopted. 
The splitting pattern of the Fe I line is a pure Zeeman triplet with 
unshifted central vr-component and two cr-components. To fit the 
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Figure 2. Measurements of radial velocity and surface magnetic field obtained for 7 Equ on 22/07/99. The Ndlll 6145.07 A pulsational RV curve (a) is 
folded with the best-fit oscillation period and is compared (c) with the variation of Bs measured from the separation of the Zeeman components of the Fe 11 
6149.25 A (centre-of-gravity line position measurements). Panels (b) and (d) illustrate the amplitude spectra for the RV and Bs (thick line - centre-of-gravity 
measurements, thin line - multiple fit of gaussians). The vertical dashed lines show the photometric pulsational periods of 7 Equ (Martinez et al. 1996). 



Table 2. Results of the analysis of RV pulsational variability of 7 Equ. The 
table gives pulsation period P, amplitude ARV and phase <f (measured 
in units of oscillation period). Superscripts "b", "r" and "c" denote respec- 
tively the blue-shifted, red-shifted and central components of the Zeeman 
resolved Fe Unes. 



Line 


P (min) 


AW (ms-l) 






22/07/99, CES/ESO 3.6-m 




Nd III 6145.07 


12.290±0.071 


454.6±23.4 


0.815±0.016 


Fell 6149.25'' 


12.290 


10.1±5.7 


0.383±0.089 


Fell6149.25"' 


12.290 


9.4±5.1 


0.215±0.085 




26/09/02, Gecko/CFHT 




Nd III 6145.07 


12.281±0.039 


165.1±10.8 


0.462±0.021 


Fell 6149.25'' 


12.281 


3.6±2.5 


0.000±0.113 


Fell 6149.25'' 


12.281 


5.2±3.2 


0.415±0.219 


Fel6173.34<= 


12.281 


5.4±14.9 


0.222±0.143 


Fel6173.34"' 


12.281 


17.2±15.8 


0.628±0.234 



relative intensity of the n- and tj-components of this line the mag- 
netic field vectors in our model have to be inclined by ~ 50° to the 
stellar surface. To look more carefully at the blending effects we 
took into account the hyperfine structure of nearby La II 6172.72 A 
and Eull 6173.03 A. An unidentified fine at A 6148.85 A, which 
shows a pulsational behaviour typical for REE lines, was synthe- 
sized with arbitrary atomic parameters to take into account its pos- 
sible blending of the blue-shifted Fe II component. 

The comparison of these calculations with the average CFHT 
spectrum is presented in Fig.Q This figure shows that the Fell 
line is free from significant blends, which enables accurate pulsa- 
tional analysis of both the red and blue components of this Zeeman 



doublet. In particular, we found that pulsational variability of the 
6148.85 A feature has negligible influence on the measurements of 
the blue component of the Fe II line. On the other hand, the blue 
cr-component of the Fe I line is blended by the Eu II 6173.03 A and 
hence is not suitable for time-resolved measurements of Bs. Con- 
sequently, RV analysis and field modulus measurements using the 
Fe I line were based on the tt- and red a-component. 

We determined the time-dependent position of the centres of 
Zeeman components using centre-of-gravity measurements within 
the spectral regions indicated in Fig.Q We note that an alterna- 
tive technique of fitting a superposition of three gaussian profiles 
to the 6148.85 A and Fell line represents a better way to derive 
Bs from time-averaged spectra because it allows to model partially 
resolved Fe II components and remove the blending contribution of 
the 6148.85 A feature. However, this method may encounter diffi- 
culties in describing time-resolved profiles of REE lines, including 
the 6148.85 A line, which show substantial RV shifts and exhibit 
profile asymmetries associated with non-radial pulsation velocity 
field. Consequently, we prefer to use the centre-of-gravity tech- 
nique throughout this paper and show the gaussian fitting results 
in Figs.|2j0for comparison purpose only. 

The RV obtained from the iron lines were compared with the 
outstanding pulsational variation seen in the Nd III 6145.07 A line. 
This feature allows us to verify the presence of rapid spectroscopic 
variability during our observations of 7 Equ and to determine the 
oscillation period and amplitude with high accuracy. For each avail- 
able dataset we established the period of the spectroscopic pul- 
sational variation using this Ndlll line and then fitted the phase 
(counted from the start HJD, see Table0 and the amplitude of the 
Fe RV measurements using a non-linear least-squares technique. 
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Figure 3. The same as Fig.l2lfor the RV and Bs determined with the time-resolved spectra of 7 Equ obtained on 26/09/02. 



Results of this analysis are summarized in Table|2| We do not find 
any pulsational variation of the RV determined for the Zeeman re- 
solved components of the Fe lines exceeding about 5-10 ms~^ for 
the Fe 11 and ~ 20 m s~^ for the Fe I line. 

The measurements of the Fell line reported in Table |2| super- 
sede tentative RV amplitude of 64 ms~^ given by Kochukhov & 
Ryabchikova ( 2001a 1, who used the same CES dataset as studied 
here. This difference mainly comes from the fact that in the present 
paper we adopt a fixed pulsation period (P — 12.29 min) for the 
analysis of the individual Fe II Zeeman components, while a set of 
four photometric periods was tested in our previous paper and the 
corresponding maximum RV amplitude for P — 11.93 min and for 
the whole line was reported. 

The measurements of the line centre positions were converted 
to surface field using the expression: 



' 4.67 X 10-"A2g' ^ ' 

where AAmag is the half of the wavelength difference between 
the red and blue Zeeman components of the Fe II line or the dis- 
tance between the red a- and the 7r-component of the Fe I line, g 
is the mean Lande factor and Ao is the laboratory wavelength of a 
line. Based on the information available from the VALD database 
(Kupka et al. 1999 1 we adopted g — 2.5 and 1.35 for the Fel 
and Fell lines respectively. The internal precision of our mag- 
netic measurements is estimated to be ~ 10 G for the Fell line and 
«70 G for the weaker Fe I line. At the same time the average field 
strengths derived from the two lines are substantially different in 
the CFHT spectra: centre-of-gravity measurements of the Fe II line 
give 3839 ± 9 G (3935 ± 14 G is obtained with multiple fit of three 
gaussians), while {Bs) = 4181 ± 66 G (4200 ± 84 G) is derived 
using the weaker Fe I line. This ~ 300 G difference is confirmed 
by the spectrum synthesis. The discrepant Bs might be related to 



Table 3. Results of the analysis of time-resolved Bs measurements for 
7 Equ. The columns give spectral lines used for Bs determination, magnetic 
pulsation amplitude ABs, phase and RMS scatter aBs of Bs - The last 
column reports an amplitude 5Bs of surface field variabihty which would 
have been detected with our data at the 3a confidence level. 



Line 


ABs (G) ^ 


(tBs (G) 


SBs (G) 




22/07/99, CES/ESO 3.6-m 






Feu 6149.25 


4.3±3.2 0.037±0.120 


13.0 


13 




26/09/02, Gecko/CFHT 






Feu 6149.25 


1.2±1.6 0.294±0.219 


9.2 


5 


Fel 6173.34 


7.9±11.8 0.628±0.234 


66.4 


28 



inaccuracy of the tabulated Lande factors, or it could be real and 
reflect different horizontal and/or vertical formation regions of the 
two diagnostic lines, or the effects of saturation in the Fe II line. An 
average field strength of 3903 ± 13 G (3941 ± 19 G) was obtained 
from the Fe II line in the ESO spectra of 7 Equ. 

Table|3|and Figs.|2j0present results of our time-series anal- 
ysis of the surface field measurements. We see absolutely no ev- 
idence for any magnetic variations during the pulsation cycle in 
7 Equ. Formal results of the least-squares fits with a fixed pulsation 
period as derived from the Nd III line indicate insignificant ampli- 
tudes, all below 10 G. Our most precise time-resolved magnetic 
measurements are derived from the CFHT dataset. Analysis of the 
6149.25 A line in these spectra results in formal Bs amplitude of 
just ~1 G, with the highest noise peaks in the amplitude spectrum 
(Fig. |3jl) not exceeding 5 G over the whole period domain typical 
for roAp pulsations. 

A 3ct upper limit of the amplitude of magnetic variability was 
estimated from Monte Carlo simulations by sampling a sinusoidal 
signal at the phases of our observations and adding a random noise, 
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Figure 4. Magnetic field measurements obtained for 7 Equ on 26/09/02 using the Fe 1 6173.34 A line. The panel (a) shows Bs phased with the best-fit pulsation 
period, the panel (b) presents the amplitude spectrum for Bs (thick line - centre-of-gravity measurements, thin line - multiple fit of gaussians). 



characteristic of the scatter of the Bs measurements. This rigor- 
ous statistical estimate is reported in Table |3| and indicates that, at 
the 3(T confidence level, no magnetic variability with the amplitude 
555 G is seen in 7 Equ during the night of our CFHT observation. 



4 CONCLUSIONS 

Our time-resolved magnetic measurements of 7 Equ have achieved 
the highest precision for a roAp star, but reveal no evidence of 
pulsational modulation of the field strength. We constrain possible 
magnetic variability to be below ~5 G in Fe lines. These results 
suggest that the marginal detection of ~100 G variability dur- 
ing the pulsation cycle reported in 7 Equ by Savanov et al. 12003 > 
is spurious, and probably stems from insufficient precision of the 
Bs measurements in that study. 

The null result reported in the present paper complements 
non-detection of the pulsational variability of Bi determined from 
Nd III lines (Kochukhov et al. 2004 (. It should be recalled that the 
Fe lines studied in our paper and the REE lines showing strong pul- 
sational RV modulation are formed at substantially different atmo- 
spheric depths due to the extreme stratification of chemical abun- 
dances in cool Ap stars and in 7 Equ in particular. Stratification 
analysis presented by Ryabchikova et al. ( 2002 ( allows us to con- 
clude that Ndlll lines sample very high atmospheric layers with 
optical depths logrsooo —7, while the Fe lines are formed be- 
low logrsooo ~ — 1. The striking difference in the RV amplitudes 
of Fe and Nd III lines is then attributed to an outward increase of 
pulsational amplitude by roughly a factor of 50-100. This increase 
is none the less not accompanied by an increase or even the pres- 
ence of any observable oscillations of the magnetic field structure. 
Combining the results of this paper with the study of Kochukhov 
et al. (20041, and taking into account that for 7 Equ Bs ~ 2.QBi, 
we find that at all observed atmospheric depths Bs changes by less 
than about 0.5 G per m s^^ of corresponding velocity oscillations. 
Hence, possible magnetic variability is constrained to be below 1 % 
of the field strength. 

We note that 7 Equ is distinguished among the roAp stars by 
its strong magnetic field, sharp lines and the exceptionally high am- 
plitude of pulsational line profile changes. Yet this star defies any 
attempts to detect magnetic variability with pulsation phase. This 
suggests that rapid magnetic modulation may be even more difficult 
to detect in other roAp pulsators. The outcome of our monitoring of 
Bi and Bs in 7 Equ demonstrates that very accurate measurements 
of the Zeeman resolved lines can yield more precise magnetic time- 



series compared with the polarimetric Bi observations. Therefore, 
the most promising direction for future attempts to detect magnetic 
oscillations in 7 Equ (which are in any case unlikely to exceed a few 
tens of gauss) would be to observe those Nd II, Nd III and Pr III lines 
which are characterized by large pulsational RV shifts and at the 
same time show Zeeman resolved profiles. Unfortunately, the extra 
broadening of the pulsating lines (see Kochukhov & Ryabchikova 
^OOlal strongly smears observed Zeeman structure and makes pro- 
posed study of Bs oscillations extremely difficult. 
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